/**********************************************************************************/
/* This class is used to get the drag force, drag acceleration, Rep for           */
/* linear motion                                                                  */
/********************************************************************************/
#pragma once
#include <cassert>
#include "Slip.h"
#include "mathParticleTracking.h"

using namespace std;
class PLinearMotion
{
public:
	PLinearMotion() { Rep = 0.0; coefDrag = 0.0; dragForce = vector<double>(2, 0.0); accDrag = vector<double>(2, 0.0); }
	PLinearMotion(const PLinearMotion& p0);
	PLinearMotion(const Slip& pSlip, const Particle& p_in, const Cphase& cphase_in);
	double Cd(double Rep);
	double Rep;
	double coefDrag;
	vector<double> dragForce, accDrag;
};

